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STEADY STATE SOLUTIONS OF TRANSMISSION 
LINE EQUATIONS 

S. O. Rice 

Methods of obtaining the steady state voltages and currents in a 
uniform transmission line consisting of several parallel wires are described 
in Part I. This line may or may not be acted upon by an externally 
impressed field distributed along its length. A square matrix r, which is a 
generalization of the propagation constant 7 for a single circuit, is intro- 
duced. Matrix expressions obtained for the voltages and currents in- 
volve r in much the same way as the corresponding single circuit expres- 
sions involve 7. In Part II similar methods are described for obtaining 
the voltages and currents in a transmission line composed of a number of 
multi-terminal symmetrical sections connected in tandem. Expressions 
for the voltages and currents in a line composed of unsymmetrical sections 
are also given. These sections may or may not contain generators. 

THE transmission lines considered here are of two kinds, namely the 
uniform transmission line, and the transmission line consisting of a 
number of identical sections connected in tandem. The problem discussed 
is that of determining the steady state electrical behavior of these lines 
when the terminal conditions are given. Often there arises the problem 
of determining the currents induced in a uniform transmission line by an 
arbitrary impressed field of some fixed frequency or of determining the cur- 
rents produced by generators placed in the branches of the sections if the 
line is of the second kind. This is the type of problem with which we shall 
be particularly concerned. 

In dealing with the uniform transmission line it is found convenient to 
introduce a matrix F, which is a generalization of the propagation constant 7 
for a single wire with ground return, or for a single circuit. This enables us 
to obtain matrix expressions for the currents and voltages which are similar 
in form to the single circuit expressions. 

A similar situation exists for the transmission line composed of a number 
of symmetrical sections. However, when the sections are unsymmetrical 
the corresponding procedure does not appear to yield a corresponding sim- 
plification and the formulas are considerably more complicated than in the 
symmetrical case. 

This paper is divided into two parts corresponding to the two kinds of 
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transmission lines. The first part discusses the uniform line. After a 
statement of the transmission equations in matrix form, expressions for the 
voltages and currents are given. Two methods of evaluating these ex- 
pressions are described. The first is based upon a property possessed by 
many transmission systems, namely that the various modes of propagation 
have nearly the same speed. The second method is based upon equations 
which may be obtained by the formal application of a theorem due to 
Sylvester. The first part concludes with the proof that these two methods 
lead to the correct results. 

After a short introduction the second part discusses the difference equa- 
tions which govern the transmission in a Hne composed of multi-terminal 
sections. The sections may contain generators. Expressions for the volt- 
ages and currents in a symmetrical section line, i.e. a line whose sections are 
symmetrical, are stated and proved in much the same order as the corre- 
sponding expressions for the uniform line. A discussion of the unsymmetri- 
cal section line concludes the second part. 

A sketch of the solution of the uniform transmission line equations by 
the classical method is given in Appendix I. In Appendices II and III 
methods are described for solving the symmetrical section line difference 
equations. These methods are similar to the one of Appendix I. The 
method of Appendix III uses section constants which may be obtained from 
measurements made at one end of a typical section. 

Part I 

Uniform Transmission Lines 

J.l Differential Eqiialions 

For the sake of convenience in writing down equations we shall assume 
that the particular line under consideration consists of three parallel wires 
with ground return, or of three parallel circuits, denoted by the subscripts 
a, b, and c respectively. The differential equations for this line in an arbi- 
trary impressed field are 

dx 

^ = -ZtJ. - Ztth - ZiJc + k(.x) (1.1) 

ax 

dx 

' These equations are given in substance by J. R. Carson and R. S. Hoyt, B.S.TJ., 
Vol. 6, pp. 495-545 {1927). Equations (1.2) are equivalent to tlieir equation (90) and 
equations (1.1) may be obtained by combining their equations (83), (84), and (94). We 
shall use the term "impressed field" to mean a field distributed along the line. According 
to our convention there is no impressed field when the line is energized only at the terminals. 
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and 



ifx 

(Hi, 
dx 

die 
dx 



= —YaaVa — Vab^b — Yac-Vc + ^oW 
= —YbaVa — YhbVb — YkoVo + tb{x) 
= —YcaVa — YcbVb — Yr<,Vc + tc{x) 



;i.2) 



where Zab ^ Z^a , I'„b = Fja , etc. If we are dealing with three parallel 
wires la{.x), lb{x), l^x) are the longitudinal components of the electric force 
of the impressed field at the wire surfaces; ta{x), tb{x), t^ix) are specified 
by the admittance of the direct leakage paths and the values of the im- 
pressed potentials at the wires. If there are no direct leakage paths the 
/'s are zero. 

In order to put these equations in matrix form^ we introduce the column 
matrices 



(1.3) 



The equations (1.1) and (1.2) may now be written as 

dv 



dx 
di 



= -Zi-\- l(x) 



5^ = -Yv + tix) 
dx 

and these are the equations to be solved. 
When there is no impressed field equations (1.5) give 



dx^ 



= ZYv 



d.v- 
2 Cf. L. A. Pipes, Phil. Mag., Vol. 24 (1937), p. 97. 
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, l{x) = 
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Jc{x)_ 
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t(x) = 


'Ux)' 

h{x) 

Mx)_ 


and the symmetrical square matrices 
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and the analogy with the one circuit case leads us to put 

p = ZF, r = Vzy (1-7) 

where T is a square matrix representing a generalization of the propagation 
constant. Putting aside for the moment the question of interpreting the 
square root, we note that interchanging the rows and columns in P = ZY 
gives 

r'2 = Y'Z' = YZ, T' = VYz (1-8) 

where the primes denote transposition. Y' and Z' are equal to Y and Z 
respectively because of their symmetry. We thus expect V to be associated 
with the propagation of i in the same way that r is associated with the 
propagation of v. 
1.2 Statement of Results for an Infinite Line—No Impressed Field 

It is shown that when there is no impressed field the voltages and currents 
at any point a: in a transmission line extending from x = to x = w are 
given by 

v{x) = e~'Mo) = e~''^ZJ(p) 

Hx) = e~'^'i(o) (1.9) 

v(x) = Zoi{x) 

wliere e^'^ is the square matrix defined by the convergent series of matrices' 

,'^r^/--^ + ^_^L+... (1.10) 

1! ^ 2! 3! ^ 

and e"""^' is the transposed of e'"'^ . I denotes the unit matrix. Z^ is a 
square matrix and is called the characteristic impedance matrix: 

z„^ v-'z^ rr-' (1.11) 

Additional expressions of the same type for Z„ are given by equations 

(1.45). The matrix e'^'^Zo, being of the nature of a transfer impedance, 

is symmetrical. 

The matrices e^^^ and Zo may be computed in several ways, the choice 

depending upon the circumstances. The first method to be described is 

useful when x is not too large and when the propagation constants of the 

various modes of propagation are nearly equal to each other. In the case 

of open-wire lines these propagation constants are grouped around the 

value yw/f where v is of the order of 180,000 miles per second. The second 

method may be used for all cases, including those for which the series in 

» Frazer, Duncan and Collar, "Elementary Matrices," Cambridge University Press, 
§2.5. In the work which follows, this text will be referred to as "F.D.C." 
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the first method converge too slowly to be of value. However, it requires 
the solution of an mth degree equation and the determination of the m 
modes of propagation where m is the number of circuits. For m = 2 
this is no handicap and the method is quite convenient. In this case 
the method is closely related to one described by John Riordan in an un- 
published memorandum. 

First Method: Multiply the matrices Z and Y together to obtain ZV. 
Choose the number 7- in 

Zr = /7- + -K, (1.12) 

where / is the unit matrix, so that the elements of R are small in comparison 
with 7^. For many transmission lines it is possible to do this. F may be 
obtained by using the binomial theorem to expand the square root in the 
formula 

r = VZV = y(I + y-^R)\ (1.13) 

where 7 is that square root of 7^ whose real and imaginary parts are non- 
negative. In carrying out the work it is convenient to introduce the matrix 
S whose elements are small in comparison with unity. 

r = 7(/ + 5) (1.14) 

To compute S, first compute the matrix R/2y^ and then use the power 

series 

^ ^ (v) " 2 (2-?) + I (2?) - I (2TV 

^ 8 \27V 16 \27V 
Tjiis series will usually converge rapidly. The matrix e'"'' is given liy 

where s is a number, z = yx, and e~^^ is to be computed from 



(1.15) 



c = c -e (1.16) 



.—IS 



zS , {%Sf (zS) 



e '^ is obtained from c~''^ by interchanging the rows and columns. The 
characteristic impedance matrix may be obtained from (1.11), 

Zo = ri-\ 

after computing V from S as in (1.14). 
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If only e~^^ is required the following series may be used. 

where R, 7, and 2 have the same meaning as above and the coefficients are 
computed from 



/ D A 

In the first term of the series ( — j denotes /. 

Second Method: V, e"^^ and Zo may be regarded as functions of the square 
matrix ZY. In order to express these functions in a form suitable for calcu- 
lation we apply Sylvester's theo^em^ The characteristic matrix of ZY is 

/(y) ^ yH - ZY {1.19} 

where now 7^ is regarded as a variable instead of a fixed number as in the 
first method. We shall suppose that ZY is a square matrix of order m 
and that the roots 71 , 72 , ■ ■ ■ 7m of the characteristic function, i.e. of the 
determinant al equation 

1/(7^)1 = 0, (1.20) 

are distinct. Let the matrix F{i^) be the adjoint of j{t) and denote the 
derivative of the characteristic function by 

i/(7^r = j^i/(TOi (1-21) 

Since y\, 72, • ■ ■ 7L are all different 1 Jiyl) j'^' is unequal to zero for r = 
1, 2, - - ■ m. Sylvester's theorem says that if P{ZY) is any polynomial 
m ZY then 

P{ZY) = i: N{yl)P{yl) (1.2Z) 



where i'(7r) is a scalar (and thus deviates from our convention that capital 
letters denote square matrices). N{y\) is a square n:atrix: 

Fiji) 

Ifiyl) 

When m = 2, N{yl) is equal to / - N{yl). 

* F.D.C. §3.9. The w and X of the reference are the ZY and 7^ of the present section. 



N{yl) = 4^]:^, (1-23) 
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Applying (1.22) to T, e~''^ and Zo even though they are not polynomials 
in ZY gives results which may be verified to be true. 

r = \/ZY = E ^(Tr)7r 

Z.= (ZY)''Y-' = j:N(yl)yrY'' ^^'^^^ 

where the summations extend from r = 1 to r — tn and 71 , 72 , ■ ■ • 7m 
are the square roots of 7? , 72 , ■ ■ ■ y„ respectively whose real parts are non- 
negative. 7i , 72 , • ■ ■ 7m are also the propagation constants of the "normal 
modes" of propagation. Some light is thrown on the physical significance 
of the matrix N{y,) by supposing that only the rth normal mode is being 
propagated on the transmission line. A'(7r) is such that it can be expressed 
as a column matrix times a row matrix. The voltages in circuits 1,2, ■ ■ ■ m 
are proportional to the first, second, ■ - ■ mth elements, respectively of the 
column matrix. The currents in circuits 1, 2, - - - m are proportional to the 
corresponding elements in the row matrix. 
1.3 Results for Any Uniform Line — No Impressed Field 

When the length of the line is finite the voltages and currents may be 
expressed as 

v(x) = cosh xT v(o) — sinh .iT Z„i(o) 

(1.25) 
i(x) = - sinh .vr' Z7'p(o) + cosh .iT' i(o) 

where Z^ and r have the same meaning as before. The matrices sinh xVZo 
and sinh xF'ZJ^ are symmetrical. The square matrices cosh xT and sinh .rr 
are defined by the series 

x'r x*r* 
cosh :vr = 7 + -^ + -^ + . . . 

smh ^r = 7-, + ^ + ■ ■ ■ 

cosh xT' is obtained by interchanging the rows and columns of cosh r^r 
and sinh xT' is obtained similarly from sinh xT. Solving {1.25} for v{o) 
and i(o) gives 

v(o) = cosh xV v{x) + sinh aT ZJ^x) 

i{o) - sinh xV Z7\{x) + cosh .xT' i{x) 

As in the case of the infinite line, we have two ways of computing the 
coefficients of v{o) and 1(0) in the expressions (1.25) for v(x) and i(x). 
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First Method: Choose a number 7^ and compute the matrices R, S, T, Zo 
as described in the first method for the infinite Hne. The matrix e' is 

given by 

xT I zS 

e = e -G 



where z = yx and e'' is computed from the series 

„2S _ J _.zS ^^z S j^ 
' -^ + l! + "2!" + "'" 

If the elements of zS are so large that the series converges slowly it may be 

(zS\ 
worthwhile to divide zS by 16, say, compute exp ( — j from the series, and 

then obtain e^ by four matrix multiplications. When e^ is known its 
inverse e~'^ can be computed and e"""^ obtained from {1.16). The hyper- 
bolic functions are given by 

coshxr=He'' + ^^'') 

• u n 1 / iT -xV\ ^ ' ^ 

smh xV = ^ {e — e ) 

which follow from the series definitions of the various matrices. 

If only the coefiicients in (1.25) are required we may choose 7- and com- 
pute R and powers of the matrix Rx/2y. Then the coefficients in (1.25) 
are given by 

cosh.r=^^^"V«.(^) 






P! 



smhxrZ.= Jli^) "^Z (1.28) 



^o\2yJ p\y 

P=o\27/ ^'7 

where R' is the transposed of R, and tlie scalar coefficient a,,(s) is a function 
of 3 = yx given by 

00(2) = cosh 2 ^1(2) = sinh s 

sinh z 



02(2) = cosh z — 



2 



(1.29) 



Qj^iiz) = Cp(z) ap+i(z), 

2 



and it is understood that {Rx/2y)° = I. 
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Second Method: Compute the propagation constants 71 , 72 , ■■■ 7n. 
and the square matrices AX7r) given by (1.23) as in the second method for 

the infinite line. Then 

cosh a:r = S iV"(7r) cosh ^7, 
sinh xT Z„ = sinh xV VY~^ 

= S N(yl) sinh ry, y,V~' (1.30) 

sinh xT' Z~^ = sinh xV V'~^Y 

= J: N'iyl) ^'"^^^- Y 

7r 

where N'{yl) is the transposed of 7^(7^), N(yl) being defined by (1.23), 
and the summations extend from r = 1 to r = ?k. 

When the transmission line consists of perfectly conducting wires strung 
on perfect insulators over a perfectly conducting earth the magnetic and 
electrostatic fields are related so as to make Z equal to yll'~^ where 

yo = jo}/c, 

u being 2x times the frequency and c the speed of light. 

It is interesting to apply the first method of solution to this line. Even 
though the proof of the first method, which is given in §1.10, does not cover 
this case there seems to be little doubt that the correct answer is obtained. 

We have 

ZY = yll 

Choosing 7 = 7« gives R = and therefore S = 0. It follows that 
r = yJ,Z,= T-' Z^yT'Z 
cosh xT = cosh (xyol) = cosh a-7o / 
sinh .-iT Zo = sinh xyoy'^^ Z 
sinh .rr' 27^ = sinh xyo yoZ~^ 

When these are put into equations (1.25) the expressions for v{x) and i{x) 
in a perfect transmission line are obtained; 

v{x) = cos)ixyov{o) — — Zi{6) 

'^0 (1.31) 

i{x) = —yosinh xyoZ~^v{o) + cosh xyoi{o) 



140 



BELL SYSTEM TECHNICAL JOURNAL 



(1.32) 



1.4 Results for Any Ufiiform Line— Impressed Field 

The differential equations to be satisfied in this case are given by (1.5). 
A solution which reduces to v{o) and tip) at .x = is 

v{x) = cosh xV v{o) — sinh xV Zoi{6) 

+ f cosh {x - Or 1(0 d^ - { sinh (x - ^TZoii^d^ 
Jo *'o 

i(x) = -sinh .vr' Z~'v{o) + cosh xT' i{o) 

- r sinh (.r - 0T'Z7'm d^ + f cosh {x - ^T'ti^d^ 

Jo Jo 

The matrices cosh xV, sinh xT and Zo are the same as the ones discussed in 
§1.2 and §1.3. The elements of the integral^ of a matrix U (U is not neces- 
sarily a square matrix) are given by the integrals of the corresponding 
elements of U. 

In many cases of practical interest the impressed field varies exponentially 
with respect to x. The column matrices l(x) and t(x) may then be ex- 
pressed as 



l(x) = e 



-xB 






t{x) = 






{1J3) 



;i.34) 



where the X's and t's are constants and 6 is the propagation constant of 
the impressed field in the direction of the line. The integrations in the 
expressions (1.32) may be performed with the result 

v{x) = cosh .vr d(o) - sinh xF Z„i(o) 

+ 1 (6^' - e-^'i) (r + on-\\ - Z,r) 
- 1 (e-^^ - fi-^V) (r - err\\ + z^r) 

i{x) = -sinh -rr' Z^^ip) + cosh xV i{o) 

+ i (/^' - e-^'l) {V + diy'w - z:'\) 
_ i (g-^^' _ ,-7) (r' - di)-\T + Z7^) 

provided that the inverse matrices exist. The matrix {e" — e ' /)(r' + 
eiy' is the transposed of (/"■ - e'''l){V + BI)~\ etc. If one of these 
matrices, say T - 61, has no inverse then it is necessary to evaluate the 

6F.D.C. §2.10. 
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corresponding integral in some other way. Thus it may be advantageous 
to use the formula 

(1.35) 
= e-'^^y:! + 1^ (r - dl) + "^ (r - dif + . . .] 

Two special cases of (1.34) are of interest. When the line is shorted at 
both ends, v{o) = v{x) = 0, where x is the line length, and 

i{o) ^ i Z-' (sinh xV)-\{c'' - e'^'mv + diy^X - Z,t) 

- {e~'^ - .-^7)(r - W)-'(X + Z,r)\ 

i(x) = ^ 27' (sinh xrrW'' - e-'lKT - BFr\\ + Z„r) 

- {€-'''- e''l){T -\- ei)-\\ - Z.t)] 

When the line is terminated in its characteristic impedance at both ends, 

v(o) = —ZJ(o),v(x) = ZJ{x),iind 

i(o) = Ul - c-''c-''"){T' + dI)-\Z-\ - r) 

i{x) = -h (£-"" - c-''l){V' - eir\Z-'K + r) 

The matrices occurring in the expressions (1.34) for v{x) and i(.v) may 
be computed by the first or second method described for the uniform line 
in the absence of an impressed field. The second method involves the use 
of expansions similar to 

{e'^ -e-''l){T-\-el)-\\ - Zot) 

(g-r -e-V)(r' + e/)"(z;"x - t) 

where the summations run from r = 1 to r ^ ?« and N'(yr) is the transposed 
of the square matrix N{yl) given by (1.23). In obtaining these expansions 
by Sylvester's theorem, Z^ in the first is replaced by F Z and Zo in the 
second by T' V. 
If we assume that an impressed field acts upon the perfect transmission 
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line of equations (1.31), we see that t{x) = because there are no direct 
leakage paths. We may also write 

(^r - e^-'DiV + dl)-' = {e'-""! - e-'' I){yJ + 0/)"' 



7.+ 
From this and similar equations it follows that 

f \ ^. r \ sinh xy^ ^ ., •. 

vix) = cosh xyoV[o) — Ityo) 

7o 



_ g-" e-^yo _ g-"- 



1 r /^° . _ . ^_ I ^ 

2L 7. + fi 7«-^ J (1.37) 

i{x) = — 7o sinh xyo Z~^v{o) + cosh xyai{o) 

yo-\-9 7o - e J 

1.5 Results for Infinite Uniform Line — Impressed Field 

When the line extends from a; = to a- = ^ and the impressed field is such 
that the voltages and currents remain finite at x = =» , the appropriate solu- 
tions may be obtained from the results of §1.4 by a limiting process. The 
condition that v{x) remain finite suggests that the coefficient of e'^ be zero 
in the expression (1.32) for v{x). This gives a relation between vip) and 
i{o) which must be satisfied: 

v{o) = ZJ{o) - [ e-^''m - Zot(0]d^ (1-38) 

If the impressed field varies exponentially with x exjiression (1.34) gives 

v{o) = ZJio) - (r + ei)-\\ - Zot) (1.39) 

Expressions for v{x) and i{x) may be obtained by using relations (1.38) 
and (1.39) in (1.32) and (1.34) respectively. As these are somewhat 
lengthy we shall state only two which follow from (1.39). 

v(x) = e~^^v(a) 

+ Ke"'' - e-"i)[(.v + eir\\ - Z„t) 

- (r - dir\\ + z„t)] ^ ^ 

(1.40) 

i(x) = ZT'e-'^'vio) 

+ ^Z-'ie-^^' + e-''l){V + dI)-\\ - Z„t) 

- ^Z7\e-' - e-^'mr - dir\-h + Z„t) 



£^^ 
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Two similar expressions may be obtained in which the initial current i{o) 
instead of v{o) appears on the right. If the line is terminated in its char- 
acteristic impedance &tx= 0, v(o) = — ZJ.(o), and the voltages and currents 
produced by the impressed field are 

v{o)= -Hr + wr(x-2.T) 

(1.41) 

i(o) = i Z7'(r + ei)-\x - Zot) 

As in §1.4 these expressions may be computed by the first and second 
methods described in §1.3. For example, the application of the second 
method to the relation (1.39) which must exist between v{o) and i(o) in an 
infinite Une gives 

v{o) = t, N{yl) \^ i{o) - -^^ (x - ^ r)1 (1.42) 

where iV(7r) is the square matrix (1.23). 

1.6 Outline of Proofs 

The proof of the results which have been stated is divided into three parts. 
In the first part it is shown that if T is a matrix such that (a) its square is 
ZY and (b) every element in the matrix e"'^ approaches zero as :ii; — » =0 , 
then the expressions for v(x) and i(x) involving F and Zo satisfy the trans- 
mission line equations. In the second part of the proof it is shown that if 
certain requirements are met T as obtained by the first method satisfies the 
conditinns (a) and (b) and hence the expressions for v(x) and i{x) given by 
the first method are correct. The third part of the proof discusses a general 
procedure which may be used to prove the equations which constitute the 
second method. 

Both the second and the third parts of the proof are based upon the solu- 
tion of the transmission line equations which is sketched in Appendix I. 
This solution assumes that the propagation constants of the normal modes 
of propagation are unequal, and our proofs are limited accordingly. How- 
ever, considerations of continuity seem to show that the first method is 
valid even when two or more propagation constants are equal. Under the 
same circumstances the second method suggests the use of the confluent 
form of Sylvester's theorem.^ 

1.7 Relations Obtained by Considering An Infinite Line 

We suppose that we are going to deal with transmission lines possessing 
the non-singular, symmetrical impedance and admittance matrices Z and Y. 
We further suppose that, by some means or other, we have determined a 
matrix V which satisfies the two conditions; (a) the square of Y is 

P = ZY, (1.43) 

and (b) every element in the matrix c~^ approaches zero as x^ oo. 

« F.D.C. §3.10. 
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Consider a line extending from a; = to a; = co , there being no impressed 
field. Viewing the line at x = as an n terminal network shows that 
there is a symmetrical matrix Zo such that v(o) = Zoi(o). Let this be 
taken as the definition of the characteristic impedance matrix Zo- We shall 
show from the diiTerential equations of the line that 

1. The voltages and currents in the infinite line are given by 

v(x) = e~'^v{o) 
i(x) = e~'^'iio) 



(1.44) 



(1.45) 



2. The matrix Zo satisfies the relations 

z^ = T~^z = zr'~^ = ri^* = i^^r' 
zr' = z~'r = r'z^' = iT^' = r'~'Y 

v{x) = ZoUx) (1.46) 

3. The matrices Zo, Z, and Y obey the commutation rules 

*(r)z„ = z„$(r') 

*(r)z = z*crO (1.47) 

F$(r) = ^{V')Y 

where $(r) is any square matrix, such as e'"'^ , representable as a 
convergent power series in T with scalar coefficients. Furthermore, 
the matrices *(r)Zo, <i>(r)Z, and F$(r) are symmetrical. 
The differential equations of the transmission line are 

^=-a. ^ =-!'». Pi=Zy' (1-48) 

dx dx dx^ 

the third following from the first two when i is eliminated. That ii(x) = 
e~''^v(o) is a solution of the third equation may be verified by direct substi- 
tution and different iation^ Since this expression for v(x) approaches zero 
as a; — » Qo and reduces to v(o) at x = 0, it represents the voltages in an 
infinite transmission line. Hence the first equation in (1.44) is true. Set- 
ting it in the first differential equation of (1.48), putting x = 0, replacing 
vio) by Zoi(o), and noting that i(o) may be regarded as an arbitrary column 
gives 

rZo = Z (1.49) 

Since T was assumed to be non-singular, Zo is equal to T~^Z. Z is sym- 
metrical and the reciprocity theorem for electrical networks requires that Zo 
^ The difEerentiation of the exponential function is discussed in F.D.C. §2.7. 
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be symmetrical, hence 



The first group of equations in (1.45) follow from this together with the 
expression T'l'"' for Z obtained from (1.43). The second group in (1.45) 
is obtained from the first group. 

The commutation rule for Zo is obtained from (1.49) together with the 
equation obtained from (1.49) by transposition. Since Z is symmetrical 

rZo = z,v', r'z, = vz^v = zj:'\ 
r"z„ = ZoV" 

and the first of equations (1.47) follow from this. The second and third 
of equations (1.47) may be obtained similarly from the relations (1.45). 
The matrix *(r)Zo is symmetrical since its transposed is Zo [*(r)]' and 
this is equal to ZJ^{T') = $(r)Zo. A similar argument applies to the 
other matrices in (1.47). 

The expression for i{x) in (1.44) may be obtained by Maclaurin's ex- 
pansion. Setting X = in the second differential equation of (1.48), 



\dx)o 



= -Yv{o) = ~YZoi(o) = -r'i(o) 



where we have used the equality between the first and last members of the 
first equation of (1.45) and where the subscript denotes the value of the 
derivative at .v = 0. Repeated differentiation gives 

ax- ax 

d^i\ „,i(di \ „,3 



and so on. Hence 



lir^'Xi^ =-''''"' 



Equation (1.46) may now be obtained by using the commutation rule 
for Zo : 

v{x) = e"'^(o) - e~^^Zoi{o) 

= Zoe~'' i{o) = Zoi{x) 
This completes the proof of equations (1.44) to (1.47). 
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1.8 Proof of Relations for Any Uniform Line — Impressed Field 

Here it is shown that if a matrix V satisfies the two conditions of §1.7 
and if Zo is the characteristic impedance matrix defined there, then the 
voltages and currents in any uniform line are given by the expressions (1.32). 
If suitable conditions are fulfilled the relation (1.38) between i»(o) and i{o) 
for an infinite line may be obtained from (1.32). 

First of all, v{x) and i{x) reduce to the required values of v{o) and i{o) 
at X = 0. All that remains to be shown is that ii{x) and Hx) as given by 
(1.32) are solutions of the transmission line equations (1.5). By substitut- 
ing (1.32) in (1.5) and using the formulas 

— cosh a;r = r sinh xT = sinh xT V 

ax , 

-;- sinh a:r = r cosh xT = cosh xT V 
ax 

which follow immediately from the series definitions (1.26) of the hyper- 
bolic functions, we obtain two matrix equations corresponding to the two 
differential equations. The terms in these equations involving via) may 
be canceled out provided 

r sinh xT = Z sinh xV Z7^ 

(1.50) 

r' cosh xV Z-' = Y cosh xV 

and these are seen to be true from (1.45) and (1.47). The terms involving 
i{o) may be canceled by a similar argument. The terms involving l{x) 
may be canceled provided 

r sinh {x - ^)r r/($) di= f Z sinh {x - ^)r'Z7'li^) d^ 

f V cosh (x - OT'Z-'m)d^ " / ^ ^°'^ ^^ ~ ^^^ ^^^^^^ 

and these are seen to be true when x in (1.50) is replaced by (x — ^). The 
terms involving t(x) may be similariy canceled. Thus we have verified 
that v{x) and iix) as given by (1.32) are solutions of the transmission line 
equation provided that the commutation rules (1.47) and the relations 
(1.45) involving Zo of §1.7 are satisfied. This is the case when T is such 
that (a) r^ is equal to ZY and also (b) every element in e"^' approaches 
zero as x —*^ «: . 

In order to establish equation (1.38) for the V of §1.7 several assumptions 
regarding the impressed field are required. Writing the hyperbolic func- 
tions in the first of equations (1.32) in exponential form and premultiplying 
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-Tx 



both sides by 2e ^ gives 



2e ' vix) = 



v(o) - ZAo) + jf e^*'[/a) - Z„/(^)]rff1 



^e-''^[vio)-^ZAo)] 



-xT 



Jo 

When X —^ ^ equation (1.38) is obtained provided that the impressed field 

and the terminal conditions at the far end are such that (a) v(x) remains 
finite, (b) the integral in (1.38) converges, and (c), the last expression on 
the right in the equation above approaches zero as a; — > <x> . 
1.9 Derivation of Equations (1.25) 

Although equations (1.25) may be obtained by setting l(x) — t(x) = 
in §1.8, it is of some interest to derive them directly. By repeated differ- 
entiation of the equations 



dv „. di 

— = -Zi, -- = -Yv 

ax dx 



(1-48) 



the second, third and higher order derivatives may be obtained. Using 
these in Maclaurin's expansion about x = gives 



v{x) = [ 



/ + ^2F + ^(ZF)^+ ... 



v{o) 



[n 



G 6 

X „„ . .V 



/ + ||ZF + |j(ZF)^ + 



...Jzi(o) 



iW = -[^ 



3 5 

X „„ , X 



/ + ^ FZ + |j {YZf + 



(1.51) 



Yv{o) 



+ 



2 i 

X ,,„ , .I- 



/ + |^FZ-|-^(FZ)' + 



...]i(. 



These series converge for all values of .r and could be used for computation 
were it not for the unfortunate fact that in most problems a great many 
terms would be required for a satisfactory answer. For the time being, 
let r be any matrix whose square is ZY. The definitions (1.26) of the 
hyperbolic functions enable us to write (1.51) as 

vi.x) ^ cosh xV v(o) — sinh :iT V~^Zi{o) 

(1.52) 
i{x) - -sinh aT' r' ^Yv{o) + cosh xV i{p) 

If in addition to being a matrix whose square is ZY, V is also such that 
every element in e'^^ approaches zero as x —^ ^c, then we may use the 
relations (1.45) for Z<, and obtain (1.25). 



.■jBiWW'aBVW rit >-j V,"'-^- 
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Incidentally, when we put ZY = ly + R in (1.51) and rearrange the 
terms so as to get a power series in R we get the series (1.28). 
1.10 Proof of the First Method 

The first method consists essentially of determining T from the series 
expansion of (1.13): 

where (-|)„ = (-i)(i)(|) •■■(«-!) when » > and (-|). ^ 1, and 
then computing Zo and the required exponential and hyperbohc functions 
of xT. From §1.7 and §1.8 it follows that the first method gives the correct 
result provided that V as determined by (1.53) satisfies the conditions: 
(a) its square is equal to ZY and (b) every element of e"^ approaches zero 
as :*: — * °o . 

These two conditions are satisfied by the matrix 

r = PGp-' (1.54) 

where P and G are matrices defined by equations (Al.l) and (A1.3) of Ap- 
pendix I, G being a diagonal matrix whose rth element is y, . For from 
(A1.9) the square of V is 

j,2 = PG^P"^ ^ ZY 
Furthermore, 

V n\ 

«! 

= PM{x)P''^ 

where M{x) is diagonal matrix (A1.5) whose rth element is e~''"'. Since 
the real part of 7, is positive and the elements of P are independent of x 
it follows that the second condition is satisfied. 

It will now be shown that PGP~^ may be expanded in the series (1.53) 

'V 

provided that 7 may be chosen so as to make all of the points f^ = — , 

r = 1, 2, . • • m, in the complex f plane lie within that loop of the lemniscate 
\f — \\ ^ \ which contains the point t = \. For then we may write 
the rth element in G as a convergent series : 

^ (-l)n(r-7^r 



(1.56) 



-7E 



n-0 n\ 
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r 2 2 
71—7 








It 







2 2 

72—7 





P~'- 


= R\ 








2 2 

73 — T _ 







and PGP~ may be written as a convergent infin ite series, the nth term 
of which contains the matrix (assuming only three circuits for the sake of 
simplicity) 



(1.57) 



where the equality follows from the definition (1.12) of R and equation 
(A1.9) of Appendix I. This series for PGP~ is exactly the same as the 
series (1.53), and this completes the proof of the first method. 

The equations (1.18) and (1.28) which are incidental to the first method, 
will now be established for the case in which the matrix F occurring in 
them is equal to PGP~ . For then we have equation (1.55) and the 
equations 



cosh xT = P 



cosh 2:71 
cosh X72 
cosh :i;73_ 



(1.58) 



sinh xT Zo = sinh xVV ' Z 
-,-1 



where sinh xT F may be expressed in the same fashion as cosh aT, the 

sinh xYr 



rth element of the diagonal matrix being 



The elements in the 



diagonal matrices occurring in these expressions may be expanded in series 



by replacing 7^ by its representation (1.56), assuming — — 1 
using 



< 1, and 



-zx/i+r __ 



' = E ? 



LI 



4/? 



K^-\{z) 



coshVl + r-E(f) ^j 
sinh z-s/X + ^ _ ■y' A^V ^ 






irZ J. / -. 



where /p_i(s) and Kp-i{z) are Bessel functions of the first and second kinds, 
respectively, for imaginary argument. Equations (1.18) and (1.28) are 
obtained when equation (1.57) and the Bessel function recurrence relations 
are used. 

* These are special cases of formulas given in "Theory of Bessel Functions," by G. N. 
Watson, page 141. 
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1.11 Proof of the Second Method 

To establish the second method we must prove the various formulas 
which are used. These formulas all involve the square matrix N{y,) de- 
fined by (1.23). 

Since N{y1) is proportional to F{yl) it follows that N{yl) may be ex- 
pressed as 

N{yl) = p.pr--- (1.59) 

where Pr is the column matrix defined in Appendix I and pr is a row matrix 
specified by pr and N{yl). Applying Sylvester's theorem to the unit 
matrix gives 

Pi 

p2 



/ = E ^(Tr) = E PrPr = Ipl , ^ , ps] 



-P3- 



where the two matrices on the extreme right are partitioned square matrices. 
From the definition of P in Appendix I it follows that 



[pi,p2,pi\ =P, 



Pi 

P2 

LPs J 



= P^' 



(1.60) 



These relations enable us to verify the equations (1.24) when T is equal 
to PGP~\ Thus for the first of equations (1.24) 



r = PGP ' = [pi,p2, pslG 



Pi 

P2 
_P3J 



= [pl , pi , p3[ 



71 Pi 

72 P2 
_73P3J 



= E ^7rP. = E ^(Tr)7r 



The second equation of (1.24) foUows likewise from the expression (1.55) 
for e-'\ 

The third equation of (1.24) follows at once from the first when we use 
(1.45), Zo = rl-^'. The fourth equation is obtained by writing 

e'^^^Z, = PMix)P~'PGP~'r-' 

= PM{x)GP~^Y~^ 

and proceeding as in the case of the first equation. 

All of the other equations connected with the second method may be 
proved in a similar way. Incidentally, the formulas obtained by the second 
method are closely related to the "special form of solution" described in 
§6.5 of F.D.C. 
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Part II 

Transmission Lines Composed of Mitlti-terminal Sections 

2.1 Introdttctory 

Some transmission systems may be regarded as consisting of a number 
of identical sections connected in tandem. The question of determining 
the steady state electrical behavior of such a system from a knowledge 
of the properties of a single section will be considered here. 

Each section will have a certain number, say ?» -|- 1, terminals on its 
left end and an equal number on the right. The case in which there are 
only two terminals (m = 1) has been completely worked out, and some 
studies of more general cases have been made. The ones which most nearly 
approach the point of view of the present paper are those due to S. Koizumi^. 

In the present work difference equations are used to solve the general 
case in much the same manner as they have been used in studying the 
two-terminal case. This approach differs from that used by Koizumi and 
throws additional light on the problem. 

In several lists of formulas, particularly in Appendix IV, I have included 
a number of results due to Koizumi for the sake of completeness. 

2.2 Transmission Equations for a Typical Section 

We consider the equations connecting the input and output currents and 
voltages for the Hth section which is shown in Fig. 1. The directions which 
are assumed for positive current tlow are indicated by arrows. The leads 
marked play a special role in that all the voltages are 



iitn) — 



Vi(n) 



in (n) 



V2(n) 



VmCn) 



— L|Cn + i) 



— i-ztnti) 



v,(n+t) 



VaCn+O 



vm(n+i) 



i-m(rn-i) 



FtG. 1 

measured with respect to them, and the currents which they carry are the 
sum of the currents flowing into or out of the remaining terminals at the 
end under consideration. In applications to transmission lines the terminals 

would correspond to the ground or the cable sheath. 
The currents and voltages shown in Fig. 1 are related by a number of 

" Archiv fiir Elect rotechnik, Vol. 33, pp. 171-188, 609-622 (1939). See also a paper 
by M. G. Malti and S. E. Warschawski, Trans. A.I.E.E., Vol. 56, pp. 153-158 (1937). 



- -i 
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sets of 2m linear equations which may be conveniently written in matrix 
form. One such set is 



(2.1) 



v(n) = Zni{n) - Zy^Kn + 1) + An) 
v(n + 1) = Z2ii(«) - Z22i{n + 1) + «"(«) 

Zn , Zi2 , Z21 , Z22 are square matrices of order m whose elements are im- 
pedances. v(n) and i(n) are the column matrices 

Vi{n) 



vUi) = 



vM 



_Wm(M)_ 



i{n) = 



ii(n) 
iiin) 

Jm(n)_ 



The column matrices v°(n) and u°(n) arise from generators which may 
be acting within the «th section. If both ends of the section are open 
circuited so that i(n) = i{n + 1) = the equations show that y(w) = 
v°{n), v{n + 1) = u^in). Consequently, if{n) and n°{n) give the open 
circuit voltages produced on the left and right ends of the »th section by 
the internal generators. If the section is a passive network then w°(«) = 
u°{n) = and they do not appear in the equations. The subscripts on 
the square matrices, the Z's, are chosen so as to preserve the analogy for 
the simple case m = 1, where the left and right ends of the section are 
denoted by the subscripts 1 and 2, respectively. 

Solving the equation (1.1) for i{n) and i{n + 1) gives 

i(n) = Vnv(:n) + Yuvin + 1) + ^(ro) 

-i(« + 1) = Y2iv{n) + Y22v{n + 1) - fin) 

where the elements of the F's are admittances and i°{n), f{n) are the 
currents produced by the internal sources when the terminals on the right 
and left are short-circuited so that v{n) = v(n + 1) = 0. 
A third set of equations is 

vin) = Av{n + 1) + Bi{n + 1) - Bf{n) 

i{n) = Cv{n 4- 1) + Di{n + 1) - Cu°{n) 

Solving these equations for vin + 1) and iin + 1) gives 

v{n +l) = D'vin) - B'H^i) + B'i°in) 

i(n + 1) ^ -Cvin) + A'ii-n) + Cv^'in) 

There are a great many relations between the square matrices appearing 
in the equations (2.1) to (2.4). These are discussed in Appendix IV. 



(2.3) 



(lA'i 



^ 
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For symmetrical sections F21 = Fjo, Y22 = Yn, Zn = Zu, Z22 = Zn 
and equations (2.1) and (2.2) become 



v(n) =^ Zni{n) — Zui(:n + 1) + v°(n) 
v(n + 1) = Zi2i(n) — Ziii(n + 1) + u°(ii) 
i(n) = Ynvin) + Y^n + 1) + i°(n) 
-i{n + I) = Ynvin) + Ynv{n + 1) - f(n) 



(2.5) 



(2.6) 



Eliminating t(") from (2.5) and v(n) from (2.6) and using, from (A4.4), 
A = ZnZ\2 = — Vvi Fu leads to the difference equations 

v{}i + 1) + v{n - 1) - 2Av{n) - B[i°(«) - f{n - 1)] (2.7) 
iin + 1) + i{n - 1) - 2A'i{n) = C[v°{n) - ii°{n - 1)] (2.8) 

Since we also have B' = B, C = C, D' = A for symmetrical sections 
equations (2.4) become 

v{n + 1) = Av{n) - Bi{n) + Bi°{n) 

(2.9) 
i{n + 1) - -Cv{n) + A'i{n) + Cv\n) 

We assume that the distribution of the sources in the branches of a sym- 
metrical network need not be symmetrical with respect to the two ends, 
even though the impedances of the branches are. 
2.3 Statement of Results for Infinite Symmetrical Section Line — Passive 
When the sections are passive the equations to be solved are, from (2.9), 

v{)i + 1) - Av{n) - Bi{n) 

(2.10) 
i(»+ 1) = -CK«) + A'i{n) 

If the tine extends from n = io n = » the solution is 

w(«) = e~" v{o) 

iin) = e~'""i{o) (2.11) 

v(n) = Zoi{n) 

where the matrix e~ is such that (a) the equation 

e"'' + e^^ = 2A (2.12) 

is satisfied, e being the inverse of c , and (b) all the elements of the matrix 
e~"^ approach zero as » -^ 00 . In dealing with sections we shall never 
have occasion to consider T itself but only its exponential and associated 
functions. The characteristic impedance matrix Z^ is defined by the rela- 
tion between the initial currents and voltages in an infinite line 

v{o) = ZJ(o) (2.13) 



v 
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A formal solution of (2.12) may be obtained by writing it as 

cosh r = ^ (2.14) 

Then 

e~^ = cosh r — sinh T 

^ A - {A^- /)* = A - {BCf 

where the square root is to be chosen so that condition (b) for e is satis- 
fied. The characteristic impedance matrix Za is given by equations (2.34) 
of which the following two are representative. 

Zo = (sinh TT'B = sinh T C"' (2.15) 

where sinh T is given by 2 sinli T = e — e 

The wide variety of sections makes it appear unlikely that there is a 
general method of determining e~^ analogous to the first method discussed 
for the uniform line. However, in some cases rapidly convergent series 
for e~^ and e'' may be obtained. For example, suppose that the elements 
of (2A)~^ are small compared to those of 2A. Then, from (2.12), 



_r 



e 



= 2A - (2A)-> - (2A)-^ - 2(2A)-s - ■ - - 



£-•' = (2A)-^ + (2A)-^ + 2(2A)-^ + • • ■ 
Again, \i A"^ - I = BC is expressible as ly^ + R where the elements of 

R are small in comparison with t^ we have (cf. equations (1.14), (1.15)) 



e' = A -\-y 






e = A —7 



/ + 



27^ 2V27V J 



Finally, it follows from a comparison of equations (2.11) and (A2.12) that 
a suitable e~^ is given by 



e-^ = PAF-\ e-'' = QAQ^' (2.16) 

where P, Q and A are the matrices designated by the same symbols in 
Appendices II and III. 

The formal application of Sylvester's theorem leads to a method of solving 
the symmetrical section line which is analogous to the second method 
discussed for the uniform line. Thus, if P{A) is any polynomial in A, then 

PiA) = E m^r)Pitr) (2.17) 
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where Pitr) is not a square matrix but a scalar and N(^r) is the square 
matrix 

iV(f.) = ^j-^,. (2,18) 

F(^) is tlie adjoint of the characteristic matrix, 

fit) = n- A (2.19) 

and fi, i'2, ■ ■ ■ f,„ are the roots, assumed to be unequal, of the character- 
istic equation 

l/(f)l = o. 

The denominator in the expression for N{^r) is the derivative of the char- 
acteristic function: 



L|l/(f) 



The formal application of Sylvester's theorem then gives 

cosh r = ^ = 2A^(fr)fr 

e-"" = A - U' - 7)^ = 2A'(r.)Xr 
e""' = S7V(f,)xr 

Z.= (sinhr)-^^ = S y i B 



(2.20) 



where A'(fr) is given by (2.18), the summations run from r = 1 to r = m, 
and Xr is related to fr through 

2f. = X, + XT' , K = ^r~ Vfl^n (2.21) 

where the sign of the square root is chosen so that [ X^ | < 1. Xr is related 
to e~^ in the same way that f^ is related to cosh r. 
2.4 Results for Any Symmetrical Section Line — Passive 

The solutions of equations (2.10) which reduce to the given values v{o), 
i(o) at n = are 

v{n) ^ cosh nT v{o) — sinh nV Zoi{o) 

_, (2.22) 

;■(«) = —sinh nV Z„ 11(0) + cosh nV i(o) 

where e~ and Z^ are the matrices of §2.3. These may be put in slightly 
different form by using the relations 

sinh nT Zo ^ Z^ sinh nV 

sinh nV Z^' = 27^ sinh nV 



•«*Aitt^.=*-. ' -v - - — '-- --^- — -^i -k!-- . -■--■-- :.o./-i'iais.*-wi.. 
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When (2.22) are interpreted by Sylvester's theorem we obtain 

v{n) = XNi^r) \l(\7" + KHo) - ^^^^ Biio)] 

(2.23) 

i{n) = SiV'(f.) [- ^4^ Cv{o) + K^r" + ^'^)iio) 

where N'{^r) is the transposed of N{tr) iind 7V(fr) is given by (2.18) and 

the summations run from r = 1 to r = m. 

2.5 Results for Any Symmetrical Section Line — Active 

When the sections contain generators the equations to be solved are 
those of (2.9). The solutions corresponding to the initial values v(o) and 
i{o) are, for n ^ i, 

v{n) = cosh wr v{o) — sinh nV ZJiip) 

+ X^ {cosh (w - p)Y Bi%p - 1) ~ sinh (n - p)Y ZoCv°{p ~ 1) 1 

""^ (2.24) 

i{n) = -sinh nT' Z;\{o) + cosh nV i(o) 

+ E {cosh (n - p)r' Cv°ip - 1) - sinh (« - p)r' Z7'5i°(^ - 1) } 

These may be simplified somewhat by replacing ZoC and ZJ B by sinh V 
and sinh T', respectively. 

The series in the above expressions may be summed when the generators 
are such that 

v\n) = e-"*i^ i^in) = e-"V (2.25) 

where fl is a scalar and i° and v" are column matrices whose elements are 
independent of n. Thus 

v{n) ^ cosh nV v{o) — sinh tiT Zoi{o) 

+ i(c"i- _ e-^n){e^ - e-H)-\Bi° - ZoCv°) 
^ i(^-"r _ e-"fl/)(e-^ - e-n)-\Bi° + Z,Cv°) 

i(n) = - sinh nV Z:\{o) + cosh nV i{o) 

+ i(e''r' _ e-"e7)(,i" - e-'I)-\Cv° - Z7'-Bi<^) 
+ ^(e-"!-' _ ,-„s/)(,-i- _ f-e/)-i(Ci)° + Z:'5t°) 

provided that the inverse matrices exist. 



(2.26) 
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We may interpret these expressions by Sylvester's theorem. For 
example, 

v{n) = Z mr) Mxr + K)v(o) - ^-^^ ^ Biio) 

where A^(fr) is given by (2.18). 

When the line extends to n = «> and the sources and end conditions 
satisfy suitable conditions we have the relation 

w 

v(o) = Zoi(o) - E e-'"'[Bi°(p - 1) - ZoCv^'ip - 1)] (2.28) 

When the impressed field is of the form (2.25) this becomes 

v(o) = ZJ{o) - (/ - e-nr\Bi° - Z.Cv") (2.29) 

provided that the inverse matrix exists. Expressions for v(n) and i{n) 
in such an infinite line may be obtained by using (2.28) or (2.29) in (2.24) 
or (2.26). 
Applying Sylvester's theorem to (2.29) gives 

K») = t A'(f ,) Q^ - -,'^. + ^-^--- V A (2,30) 

r-l \Xr — Xr Xr — 6 2(Xr — 6 ) / 

The last term within the braces may be replaced by 

2BCv° 

(xr' - X,)(X7^"- f^) 

2.6 Derivation of the Properties of an Infinite Line 

We shall consider a symmetrical section line which is specified by the 
equations 

v{n + 1) = Av{:n) - Bi{n) 

(2.10) 
i{n + 1) = -Cv{n) + A'i{n) 

From these equations and the relations A- — BC = I, AB = BA', A'C = 
CA of (A4.6) it follows that 

v(n + 1) + v(n - 1) = 2Av(n) 

(2.31) 
i{n + 1) + i{n - 1) - 2A'i{n) 






158 BELL SYSTEM TECHNICAL JOURNAL 

If e~^ is a matrix satisfying the conditions of §2.3, namely, (a) e~^ satisfies 
the equation 

2 cosh r ^ / + e"'' = 2/1 (2.14) 

and (b), every element in e""^ approaches zero as h — * » , and if Z^ is 
defined by v(o) and i{o) for an infinite line as in (2.13), then 

1. In an infinite line 

v(n) = e~"^v{o), 

_ , (2-32) 

i(«) = e " i((j), 

v(:n) = Z,i{n) (2.33) 

2. The characteristic impedance matrix Z^ is given by 

Z„ = (sinh ry'B = ^(sinh r')^' - C"' sinh r' = sinh r C-' 
Z;' = B-' sinh r = sinh Y' B'' = (sinh T'Y'C = C(sinh T)"' 

3. The matrices Zo , -5 and C obey the commutation rules 

*(e')Z, ^ Z,*(/') 

$(/)_B = B^{e'") (2.35) 

C$(/) = $(e'")C 

where $(c^) is a square matrix representable as a sum of powers of e . 
The matrices ^(e^) Zo, *(e'') B, and C$(e^) are symmetrical. 

To prove these statements we proceed as follows: By direct substitution 
into (2.31) it is seen that v{n) = c'^^vio) is a solution by virtue of condi- 
tion (a) satisfied by e"^. Since, by condition {h),v{n)^oas «-> « it fol- 
lows that v(n) is the voltage in an infinite line. Similarly, i{7i) = e^ i{o) 
is the current in such a line. Substituting the expressions (2.32) for i)(«) 
and i{n) into the difference equations (2.10), setting n = 0, using the 
definition of Z^ , and regarding v(o) and i{o) as arbitrary columns gives 

.-^ = A - Bz;' 

(2.36) 



-r' 



e 



= -CZ, + A' 



Applying condition (a) in the form of (2.14) to these equations gives 

BZ7' = sinh r, CZ, = sinh V (2.37) 

Since the sections are symmetrical, B and C are symmetrical matrices, and 
from the reciprocal theorem for networks it follows that Zo is also sym- 
metrical. These remarks and (2.37) lead to (2.34). Setting the expres- 
sions (2.32) for vin) and i(«) in the second of the difference equations (2.10) 
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using the definition of Z^, , and regarding i{o) as an arbitrary column gives 

(A' - e-'V'' - Cc-^Z, 

Replacing A' — e~^ by CZ^, as follows from the case n = 0, and pre- 

multiplying by C~^ gives 

and this leads to the first of equations (2.35). From (2.34) and the rela- 
tions AB = BA', CA = A'C we have 

sinh T B = B sinh V cosh V B = B cosh V 

C sinh r = sinh V C C cosh V = cosh r' C 

Addition and subtraction leads to 

from which the last two of equations (2.35) follow. Since each of equations 
(2.35) expresses the equality of a matrix and its transposed, it follows that 
the matrices are symmetrical. 

Equation (2.33), which is almost self-evident on physical grounds, fol- 
lows from 

v{,i} = 6""^") = e-^^ZJio) 

= Zoc"" i{o) = Zoi{n). 

2.7 Proof of Relations for A ny Symmetrical Section Line 

The expressions (2.24) for v{n) and /(») in a hne whose sections contain 
generators may be verified to satisfy the difference equations (2.9). The 
expressions (2.34) for Zo and the commutation rules (2.35) for B and C 
are used in the verification. Setting n ^ 1 in the expressions for v{n) and 
iiji) gives the difference equations (2.9) and lience v{n) and i{n) are the 
solutions whidi correspond to the initial values v{o) and i{o). 

In order to derive the relation (2.28) between lio) and i{o) for an infinite 
line we put the hyperbolic functions in the expression (2.24) for v{n) in 
exponential form and multiply through by 2c~"^ 

2c-"^■(H) = v{o) - ZJ{o) + i: e-''^[Bi°{p - 1) - ZXv'^ip - 1)] 

+ e-"'Mo) + ZJio)] 

+ .-"' E e-'"-''^[Bi%p - 1) + ZXv%p - 1)] 

p-1 



" \ -^--tJ ^ i *^ faa **'J ■ - JAr*-'' ■*^-- 
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Hence, letting « -^ « and using condition (b) satisfied by e , equation 
(2.28) is obtained provided that (i) the terminal conditions at the far end 
are such that vin) remains finite, (ii) the sum in (2.28) converges, and (iii) 
the expression in the last line in the equation just above approaches zero 

as « — » «> . 

The results obtained by the formal application of Sylvester's theorem 
may be verified by using the results of Appendix II and writing Ni^r) as 
the product of a column matrix and a row matrix. They may also be 
verified more directly. For example, setting 7i = Oin the expressions (2.23) 
for v{n) and i{n) in any passive symmetrical section line and using 

m 

E A^(r.) = ^ (2.38) 

which follows from Sylvester's theorem, we see that v{n) and i{n) reduce 
to the appropriate values v{o) and i{o) at n = 0. Substituting v{n) and 
i{n) into the difference equations (2.10) and using 

BC = A^ - I 

in. - A)N{tr) = N{i:r)(nr ' A) = ^^ ^^^ 

shows that they are solutions. The second of the relations (2.39) follows 
from the fact that A^(fr) is proportional to the adjoint F{tr) of /(fr). In 
the third and fourth relations 

which is in agreement with the definition (2.18) of Ni^r). To establish 
the third relation we start from, " 

(tr-A)F{t) ^/|/(r)l 

(f/-^)A^(f) = /|/(f)l/l/(f)l'" 
Postmultiplication by B gives 

{^i-A)N{t)B^B\M)\/\mr 



We also have 



'"F.D.C. §3.5. 



(ti-A')F'i^) = i\m\ 
{n-A')N'i^) = i\m\/\mf' 
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Premultiplication by B and use of BA' = AB gives 

(^i-A)BN'{o = B[m\/\mr 

Hence, the third equation in (2.39) holds except possibly for i" = Tr , and 
from the concept of continuity it holds there also. The fourth equation in 
(2.39) may be proved in the same manner. 

The expression (2.20) for Zg may be obtained by letting n become very 
large in the expression (2.23) for v(n). v{o) and i(o) must be related so that 
v{n) remains finite. Since | Xr | < 1 and the X/s are unequal the coefficients 
of X7" must vanish. This requires 



\7 - X*- 

Summing r from 1 to m and using (2.38) gives the required expression for Zo . 
2.8 The Unsymmeirical Section Line 

The method used here is analogous to those described in Appendices I and 
II for the uniform line and the symmetrical section Line. The other methods 
apparently do not lead to the simplification which occurs in the symmetrical 
case. 

Equations (2.2) and (2.1) lead to the difference equations 

Ynv{n + 2) + [Yn + Y^iHn + 1) + Y2iv{n) = -i\n + 1) + f{n) (2.40) 

Zxiiin + 2) - [Zn + Zn]i{n + 1) + Zni{n) = v°{n + 1) - «''(«.) (2.41) 

Both of these equations are of the form 

Gx(n + 2) + Hxin + 1) + G'x{n) = g(«) (2.42) 

in which G and H are square matrices of order m, H being symmetrical and 
G' being the transposed of G. When the sections are passive equations 
(2.40) and (2.41) become 

Yuv(n + 2) + [I'n + Y^M^ + 1) + r2iv(n) = (2.43) 

Zui{n + 2) - [Zn + Z22]i{n + 1) + Z^iiin) = (2.44) 

In the passive, unsymmetrical case the expressions for v{?i) and i(n) are 
of the form 

v(7i) = PA"a + PA-"d 

(2.45) 
i(n) = QA"a — QA-"d 

Comparison with (A2.8) shows that in the symmetrical case P = P and 
Q ^ Q. The minus signs over P, Q, and a indicate that they are associated 
with propagation in the negative direction. The propagation constants of 
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the m modes of propagation are the same in the positive as in the negative 
direction, as indicated by the appearance of A" and A"" in (2.45). Cor- 
responding to any given propagation constant say Xr , there are two modes of 
propagation, one in a positive direction and the other in the negative 
direction. The distribution of the voltages corresponding to these two 
modes are given by the rth columns in P and P, respectively. The fact 
that P and P differ shows that the distributions difier according to the 
direction of propagation even though 'the propagation constant is the same. 
A is still the diagonal matrix defined in (A2.3) but now the computation of 
the elements Xr is more difhcult than when the section is symmetrical. 
They are defined as the roots of the equation 

\Gy+ H\-\-G'\ = (2.46) 

which are less than unity in absolute value. The second of the equations 
(A4.5) shows that the roots of (2.46) are the same whether the Z's or the 
F's are used in place of G and H. Of course, this is to be expected on 
physical grounds. The third of the equations (A4.5) may be used to show 
that the roots of (2.46) are also the roots of 

= (2.47) 

\C \D - I 

From the form of (2.46) it follows that if K is a root so is ^7^ This fact 
may be used to simplify the determination of X, . When the substitution 

2t ^ X + \-\ X - f - >/F^^l (A2.4) 

is made equation (2.46) may be written as 

= 1(G + G')t + // + (C - G) VF^n I 
= 1 (G + G')r + H I 

_]_ (j-'^ _ 1) times the sum of — -^y — determinants each ob- 
tained by replacing two columns of |(G + G')r + H\hy the cor- 
responding columns of (G — G') 

, . , r '«(w — l)(m — 2)(m — 3) , , _. 

+ (f- - 1)- times the sum of -^ ^^-j^ — determi- 
nants each obtained by replacing four columns of |(f7 + G') f + 
// |by the corresponding columns of (G — G') 

+ ... 

Tlie last equation is a polynomial of degree m in t which is to be solved 
for its roots Ti , f2 , ■ ■ ■ i"m ■ For simplicity we assume that these roots are 
distinct. X, is then determined from f^ by the relations (A2.4), the sign 
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of the radical being chosen so that | \r | < 1 as in the symmetrical case. In 
his second paper Koizumi has given a procedure which amounts to an 
alternative method of determining A. 

We shall first assume that the T's are known and that our equations are 

i{n) = Ynvifi) + Yuviti + 1) 

(2.48) 

~i{n + 1) ^ I'2iy(«) + l'*22t'{« + 1) 

As described above A may be computed from tlie determinantal equation 

|/(X)|= 

where f{\) represents the matrix 

/(\) = Vn\' + (Fu + F2o)X + Vn (2.49) 

Let Pr be proportional to any non-zero column in F{\r) where F(\) is the 
adjoint of /(A) and let Pr be proportional to any non-zero row of F(Xr). 
Then the matrices P and P in the expressions (2.45) for v(n) and i(w) are 
given by 

(2.50) 

where pr is the column obtained by transposing the row p'^ . The matrices 
Q and Q are obtained from F and P by means of the equations 

Q = VuP + rnPA - -V22P - V21PA-' 

- - , - _ (2.51) 

Q - -VnP - 1',2/^A-' ^ J\,P + FaiPA 

which are derived from (2.45) and (2.48). 
The properties of the individual columns of P and P lead to the relations 

Y12PA' + (I'u + y22)PA + y2iP - 

- . _ , _ (2.52) 

Fi^A - + (I'll + ]'-)PA-' + l\iP = 

and these guarantee that the difference equations (2.48) will be satisfied 
when the expressions (2.45) for v(n) and i{n) are used. 

When the Z's are known instead of the F's the procedure is much the 
same. The difference equations are 

v{n) = Zni(ii) — Zi2t(« + 1) 

(2.53) 
vin + 1) = Z2,i(«) - Z22i(n + 1) 

and the equation to determine the X/s is 

i/(X)l =0 



,tfj>M^'- 



. JlA»«-^< 



-1 

(2.56) 



164 BELL SYSTEM TECHNICAL JOURNAL 

where now /(X) represents the matrix 

/(X) = Z,2X= - (Zn + Z22)\ + Zn (2.54) 

Let qr be proportional to any non-zero column in F(Xr) where F(\) is the 
adjoint of /(X) and let gl be proportional to any non-zero row of F{\r). 
The matrices Q and Q in the expressions (2.45) for v(n) and j(») are given by 

where qr is the column obtained by transposing the row Qr . From (2.45) 
and (2.53) equations for P and P in terms of Q and Q are obtained: 

i' = ZnQ - Z12QA = -Z22Q + Z21QA 

P = -ZnQ + Z12QA"' = ZsaQ - Z21QA 

The difference equations (2.53) are satisfied by our expressions for v{n) and 
»(») by virtue of the relations 

ZnQA' - (Zn + ^22)^^ + Z21Q =0 

ZuQA^' - (Zn + Z22)0A"' + ZnQ = 

which are a consequence of the properties of the individual columns of 
() and Q. 

If the system extends to m = + ^ and if the voltages and currents are to 
remain finite at m = =» the elements of a must be zero and the expressions 
(2.45) for v{v) and i{n) become - 

v{7i) = PA^a = PA^P'^vio) 

i{n) = QA"a = QA''Q-^Hio) (2.58) 

v(n) = PQ~'i{n), i{n) = QP~\{n) 

where we have assumed that P~^ and Q~^ exist. We accordingly introduce 
the characteristic impedance and admittance matrices Z^ and l'„ associated 
with propagation in the positive direction, i.e., in the direction of n in- 
creasing. 

v{n) = Z„i{n), i{ii) = Yovin), Z^ = F7 
Z„ = PQ:' = Z„ - ZnQAQ-' = -Z22 + ZnCA"'(?"' (2.59) 
F„ = QP-' = I'll + YuPAP-' = -Y22 - Y2iPA-'p-' 
Incidentally, since Z„ must be a symmetrical matrix the above equations 
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show that ZuQAQ~^ and ZziQAT^Q'^ are symmetrical. Zo and Yo satisfy the 
relations 

ZjCZ„ -j- ZJ)- AZo~ B = y^Y^ ^y^ - j)Y^^ c =0 

(Zii + Z„)^'(Z„ - Z„) = Z21 , (Fi.2 + Yo)YT2{Yn ~ Y,) = Fn (2.60) 

Z^QKQ'' = P\F~''Zo YJ'AF^' = QAQ'^Y^ 

The characteristic and admittance matrices ^o and Yo associated with 
propagation in the negative direction are introduced in a similar way. 
Suppose the system extends to m = — =» . Then a = and 

i(h) = QA.-"'d = -QA-"Q~'j(o) (2.61) 

v(n) = -PQ~\{n), i{n) = -Qp-'vin) 
Hence we write 

v(n) = —Zoi(n), i{n) = —YoV{n) 

2o = PQ-' - -Zn + Zi2QA~'Q-' ^ Z22 - Z21QAQ-' (2.62) 
y = Qp-' = -^ Fn - YnPA-'p-' = Y^ + FaiPAP"' 

Zo and Yo saiisfy the relations 

ZaCZ, - ZJJ + AZ^ - B = y,BY, - Y,A + D?^ - C = 

(Zu - Zo)Z^i\Z22 - Zo) = Z12 (Kn + F„)Fri'(F22 - f „) = F12 



(2.63) 



The fact that Q'{Zo + Z„)Q = P'(F„ + F„)P is a diagonal matrix may 
be used as a check on computations. 

When the expressions (2.45) for v{n) and i{n) are placed in {2.3), f{n) and 
«**(«) being zero, we obtain the relations 



, -1 



PA ' = AP + BQ PA - AP - BQ 

(2.64) 
QA^ = CP + IX> QA = ~CP + DQ 

When the typical section contains generators the difference equation to 
be solved is of the form (2.42) 

Gx(n + 2) + IIx{n + 1) + G'x(n) = g(n) (2.42) 

This is true for symmetrical as well as unsymmetrical sections, G being a 
symmetrical matrix in the former case so that G' = G. The expressions for 
v{7i) and i{n) are those of (2.45) with the particular solutions added: 

v{n) = PA"a + PA"''a + u(n) 

(2.65) 
i{n) = QA^'a - QA'^d + j(n) 
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where P, P, Q, Q are determined as before and m(m) and j(n) depend upon 
the generators. 

Here we shall consider only the physically important case in^which the 
voltages of the generators in the nth section are proportional to e " where 6 
is a constant. In this case g(n) may be expressed as 

g{n) = ge-"' (2.66) 

where g is a column matrix whose elements are independent of n. A 
particular solution is obtained by assuming 

x{n) = ye 

Setting this in (2.42) gives 

{Ge-'' + He-' + G')y = g 

Hence a particular solution is 

x{n) = {Of-'' + He-' + G'The^' (2.67) 

This method fails when d is equal to one of the roots Xi , - ■ ■ Xm , Xi , • • ■ Xt™ . 
In this case, a particular integral may be obtained by a method similar to 
one described in §5.11 of F.D.C. 

APPENDIX I 
Classical Solution of Uniform Transmission Line Equations 

For the sake of convenience we again assume that there are three circuits 
in the transmission line. The equations to be solved are: 

^= -Zi, f=-Yi (1.48) 

dx ax 

We adopt here the notation associated with equations (1.19) and (1.20), 
f{y) being the characteristic matrix of ZY, F{y) the adjoint oi f{y), and 
7i , 72 , yli-^ ^ 3) being the roots, supposed distinct, of \f{y^) \ = 0. The 
propagation constants 7i , 72 , 73 are those square roots of 71 , 72 > 73 which in 
physical systems have a positive real part. 

The solution may be constructed" as follows: Let the column pr be 
proportional (with any convenient constant of proportionality) to any non- 
zero column of 7^(7^). The non-zero columns of ^(7^) are proportional to 
each other according to a theorem in matrix algebra. ^^ Construct the 
square matrix P from the columns pi , p2 , Pi '■ 

P= {pi,Pi.Pz\ (Al-1) 

" The method is that described in F.D.C. §5.7(a) and §5.10 
" F.D.C. §3.5 Theorem D. 
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and obtain the square matrix Q from P: 

Q - Z~'PG = YPG-' 

where G is the diagonal matrix 

~7i 
G = 72 
_0 73_ 

The voltages and currents at any point x are 

v(x) = PM(x)a + PM(-x)d 
i(x) = QM{x)a - QM{-x)a 
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(A1.2) 



(A1.3) 



(AI.4) 



where a and a. are arbitrary column matrices associated with propagation 
in the positive and negative directions of .v and M{x) is the diagonal matrix 



M{x) = 



-711 











-y,. 











g-.,* 



CA1.5) 



The values of a and a are to be determined from the boundary conditions. 
When the line extends to .t = « 



v{x) = PM{x)P~\{o) ^ Z,i{x) 
i{x) = QM{x)Q~H{o) 
where the characteristic impedance matrix Zo is given by 
Zo = PQr^ = PG^^P'^Z = PGP~^Y~^ 
= ZQG~^Q~' - y'QGQ'^ 



(A1.6) 



(A1.7) 



Since v = PrC*'' and i = q^^'' , where qr is the rth column of Q, are solutions 
the differential equations give 

{l-il - ZV)Pr - 0, (lyl - YZ)q, = 



and from these it follows that 



P'^ZYP = Q'^YZQ = G' 



(A1.8) 



(A1.9) 



The relations (A1.8) may be used to prove the following: 

1. The elements in the rth column of Q are proportional to those in the 
non-zero rows of ^(7r)- 

2. The matrix P'Q is a diagonal matrix and from this it follows that if ^ 
is any diagonal matrix 

WP->)' = Q^Q'' (ALIO) 



- - ■gitirWf i-^^^^'^ 
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3. The characteristic impedance matrix Zo satisfies the relation 

Z = Z^YZs, (Al.ll) 

4. The inverse matrices P~^ and Q~' always exist if 71 > 72 , 73 are distinct. 

APPENDIX II 

Classical Solution of Symmetrical Section Line Equations — I 

The method of this section is very similar to that of Appendix I. The 
equations to be solved are (2.10) or one of the sets 

v{n) = Zni{n) - Z^iin + 1) 

v(n -f- 1) - Zi2i{n) - Zni(n -\- 1) 

i{n.) ^ Ynv(n) + Ynvin + 1) 

-i{n + 1) - Yuv{n) + Fiij;(« + 1) 

which are obtained from (2.5) and (2.6). We shall use the notation asso- 
ciated with equation (2.19), /(f) being the characteristic matrix of A, f (f) 
the adjoint of /(f), and f 1 , fa , ■ ■ ■ fm the roots, assumed unequal, of the 
characteristic equation |/(f) ] - 0. The diagonal matrices A and 2 are 
defined by 



(A2.1) 



(A2.2) 



A - 



"Xi 
M 





X™ 



2 = 







Vf?-i •• ■ 

Vfl - 1 

VfS. - 1- 



(A2.3) 



where 



2f. - X. + \7\ Xr = fr - Vr? - 1 = 



(A2.4) 



f . + Vf? - 1 

In general, electrical energy will be dissipated in the typical section and 
from the physical significance of Xr , as seen from equations (A2.8) below, it 
follows that the sign of the radical in (A2.4) may be chosen so that | Xr | < 1. 
Since VtT^^l = tr - K = MX7' - X.) it follows that 

S = KA"' - A) (A2.5) 

Let the column matrix Sr be proportional to any non-zero column in 
F{^r) where i^(f) is the adjoint of /(f). (It follows from the theory of 
matrices that the non-zero columns of F(f r) differ from each other only by a 
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multiplying factor.) The matrix S is then formed by taking .Ji to be the 
first column, Sa the second and so on. 

S =[si,$2,--- sj (A2.6) 

Similarly let the row matrix (.^ be proportional to any nonvanishing row of 
F{^r) and form the matrix T where 

T =[h,l2,---i^] (A2.7) 

in which ir is the column matrix obtained by transposing tl ." 

Solving our difference equations for the passive case by the customary 
method gives the expressions 

v(ii) = PA"a + PA~"a 

(A2.8) 
i{n) = QA"a - QA "a 

for the voltages and the currents. P and Q are square matrices and a and a 
are column matrices whose elements are determined by the boundary con- 
ditions, a and a are of the same nature as constants of integration. The 
minus sign over a indicates that it is associated with propagation in the 
negative direction, i.e., in the direction of n decreasing. 

P and Q may be chosen in a number of ways, each choice requiring 
different values of a and a to represent the same system. In all cases, 
however, the fth column of F may be expressed as arSr where a, is a scalar 
multiplier which may depend upon r. Similarly the rth column of Q may be 
expressed as /?,/, . When either P or Q has been chosen the other one is 
fixed since equations (A2.2) and (A2.1) require 



Q = YuP + VizPA = -YnP - VnPA^' 

P = ZuQ - ZuQA = -Zi,Q + ZnQA~' 
Some useful choices are, 

1. p = s, Q ^ -YuSz = jr^s::: 

2. P = SX Q= ZV2S = CS 

3. Q = T, P = Z12TX = C'^TX 

4. (2 = rs P = -YuT - BT 



(A2.y) 



(A2.10) 



The particular choice to be made depends upon the system of difference 
equations which is being used. In choices 1 and 2, T is not required and in 
3 and 4, S is not required. However, if both S and T are known some of 

"Methods of determining Sr and tj. are available. A description will be found in 
F.D.C. §4.12. 
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the matrix multiplication may be avoided. Taking choice 1 as an example, 
we may determine the rth row of Q from the expression ^Ar ■ To determine 
pr only one element in the rth column of — ri252 need be known, for ^r is 
the quotient obtained by dividing this element by the corresponding element 
in tr . The product P'Q must be a diagonal matrix, and the same is true of 
S'T. This may serve to check computations. 

That the expressions for v(n) and i{n) given by (A2.8) and (A2.10) satisfy 
the transmission equations (A2.1), (A2.6) and (2.10) may be verified by 
direct substitution and use of 

S(A -I- A"') ^ 2.45 T{A + A'^ = 2A'T (A2.U) 

These relations follow from the properties of the individual columns of 
5 and T. 

Wlien the system extends to w = =00 must be zero in order that the 
voltages and currents may remain finite. This is true because Xr is chosen 
so that I Xr I < 1. From equations (A2.8) it follows that 

v{n) = PA"fl = PAT'^vio) 

i{n) = QA.^'a = QA^Q'^Ho) (A2.12) 

v(n) = PQ-H{n) i(n) = Qp-'v(n) 

the reciprocal matrices Q~'' and P~^ always exist when the sections are 
symmetrical and the roots f 1 , f2 , ■ ■ ■ fm distinct. The last equations in 
(A2.12) suggest the introduction of the characteristic impedance and 
admittance matrices Zo and Yo '■ 

v(n) = Zoi{n), i(n) - Y^v{n), Z„ - I^^ . 

Zo = PQ-' = Zn - ZuQAQ-' = -Zn + ZiaQA-'Q-' 

= 5S"'5"'B = S-LS~'Zri 

= ZnTZT'' = BTir'T~' (A2.13) 

F„ = QP-' = Fn + Yv.PAP-' = - Yn - Y,,PA-'p-' 

= - YuS^S'^ = CSX~'S~' 

= T^~^T~^C - -T^T~^Yn 
Not all of the expressions for Zo and Yo obtainable from (A2.10) have been 
included in {A2.13). Z^ and Y„ are symmetrical matrices. Although P and 
Q are arbitrary to some extent the same is not true of Zo and Yo . Com- 
puted values of Zo and Y^ may be checked by use of the relations 

A'-I^ (Z^ZT^f = (YiiY;)' 
ZJZZo = B, Y,BYo = C (A2.14) 

YoZu= -YnZ^ 
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Sometimes it is desirable to terminate a line consisting of a finite number 
of sections by a network which simulates an infinite line. As is known, the 
elements in one such network may be obtained from I'o . Every terminal 
is joined to every other terminal, including the return terminal (denoted 
by the subscript o), by the branches of this network. The admittance of 
the branch connecting terminal ; to terminal y, i ^ 0,j 9^ 0, is —ya where 
yij is the element in the ith row and jth column of Vo ■ The admittance 
of the branch connecting terminal i to terminal o is yn + y.2 + ■ ■ ■ + yu 
+ ■ ■ ■ + Jim , i.e., it is the sum of all the elements whose first subscript is i. 

APPENDIX III 
Classical Solution or Symmetrical Section Line Equations— II 

When the electrical properties of a typical symmetrical section are to be 
determined by measurement, equations (A2.1) and (A2.2) show that Zu 
and Fii may be obtained by measurements at one end. In order to obtain 
Vu and Zi2 measurements have to be made at both ends. Expressions for 
v(n) and i(n) will now he given which depend only upon Zn and Yn and 
hence are useful in case the measurements are restricted to one end. 

The method is based upon the equations 



v(n + 2) + v{n) = Zn[i{n) - ^n + 2)] 
i(n + 2) + i(n) = Vn[v(n) - v(n + 2)] 



(A3.1) 



which may be derived from (A2.1) and (A2.2). Combining these equations 

leads to 

[/ - ZnVuMn + 2) + v(n ~ 2)] + 2 [/ + ZnViMn) = 
[/ - YiiZuMn + 2) + i{n - 2)] + 2 [/ + YnZn]i{n) ^ 
The first step in the solution is to solve the equation 

1 m/ - ZnYn 1 - (A3.2) 

for its roots mi , M2 , ■ ■ ■ Mm which we shall suppose are distinct. The diag- 
onal matrices M and M^ are defined by 

111 ... 

M2 



M = 







|Um_ 



^,\ 



M* = 





i 
JU2 











(A3.3) 



and A is defined as in (A2.3) where X, is given bj 



y-r = 



1 + x 



1 - X 



;i 



(A3.4) 



':.iiiaii€e-ima£i<S'iS^i^i^-!^ ^ 
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The sign of fil is chosen so hat | Xr | < 1, and this is the value to be used 
in K*. However, there is -n ambiguity in the sign of Xr which is inherent 
in this method. A relation between A and M is 

M^ = (1 + A')(I - A'y' (A3.5) 

Let Ur be proportional to any non-zero column and w'r be proportional 
to any non-zero row of the matrix adjoint to [nJ — ZnYn] and form the 
matrices 

U = [wi , «2 , • ■ ■ u„,] 

W = [Wl ,W2j • • ■ Wm] 

(cf. equations (A2.6) and (A2.7) for 5 and T) where w,. is the column 
obtained from Wr- 
The voltages and currents are given, as before, by 

vin) = PA^'a + PA~"d 
^ ' _ {A2.8) 

i{n) = QA"a - QA "a 

and there is again a number of ways in which P and Q may be chosen. In 

all cases the rth column of P may be expressed as Oriir and (he rth column of 

Q as 0rWr . The equations fixing Q when P is chosen and vice versa are, 

from equations (A3.1) 

n = YnPM~^ 

(A3.6) 

P = ZnQM~' 

where -M~* is the inverse of f.f^. Equations (A3.6) may also be obtained 
from (A2.10). 

Suitable choices for P and Q are 

1. P= U, Q= YnUM^^ = ZTiUM^ 

2. P = UM\ Q = VnU = ZTiUM 

3. Q=W, P= ZnWM^^ = Yn'WM^ 

4. g = wm\ p = ZnW - ytIwm 

P'Q and U'W must be diagonal matrices. That the expressions for y(K) 
and i{n) just derived satisfy the difference equations (A3.1) may be verified 
by making use of 

UM - ZnYiiU, WM = YnZnW (A3.8) 

Equations (A3.8) follow from the properties of the individual columns of 
U and W. The characteristic impedance and admittance matrices are 



(A3.7) 
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given by 

Z^ = PQ-' = ZnQM'^Q-' = PM^p-'\Ti 

= UM-^ir'Zn = uu^ir'YTi 

= ZuWAT^W^' = YTiWM^W-"^ 
F„ = QP^' = YnPMT^P'' = QM^Q-'zr' 

= ViiUM-^ir' = zTiUM^ir' 

= W^iT^H-^'r,! = WM^W'^ZTi 

The matrices Zo and Y„ may be checked by means of the relations 

ZoVn = ZnV, , Z.YuZ, = Zxi , Y^ZuY, = r„ (A3.10) 

Another set of solutions may be based upon the equations 

2v(n) = -ZiS(n+ 1) - i{n - 1)] 

(A3.11) 
2i(n) - YuHn + 1) - v{n - 1)] 

which are derivable from {A2.1) and (A2.2}. Combining these equations 
gives, upon using YT2ZT2 = —BC, 

v{n + 2) - 2v{n) + i;(« - 2) = ^BCv{n) 

(A3.12) 
i{n + 2) - 2i{n) + i{n - 2) = 4CBi(H) 

However, we shall not consider these equations here beyond pointing out 
that they lead to 

P = Zx,Q^, Q = -Fiji's 

PS' - BCP, gs' = CBQ 

which may also be derived from (A2.10). 

APPENDIX IV 

Relations Between the Square Matrices of 
A Multi-Terminal Section 

When the reciprocal theorems of network theory are applied to equations 

(2.1) and (2.2) it is found that Zn , Z2- , I'n , I'aa are symmetrical and 

Z21 = Z'n , Y21 = Y[2 (A4.1) 

i.e., Zn and Y21 are the transposed matrices of Zn and Yu , respectively. 
Solving equations (2.1) for the currents and comparing the result with 

(2.2) shows that 

\Zn Z^y ^ \Yn F1.I f ^(n) 1 _ JYu Y^^lVv^in)! 

IZn Z22J lYn F22J' l~r{n)j lYn Y,,]lu%n)] 



I 

: — l: j_. Ill-- --■i 



^22JLF21 F22J LO ^_ 



{A4.2) 
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These are partitioned matrices. The square matrices have 2m rows and 
columns and the column matrices have 2m elements. The first of these 
relations may be written as 
Zii Z\ 

_Z21 Z2 

where / denotes the unit diagonal matrix of order m. Multiplying the two 
matrices on the left together and equating the elements of the product 
to the elements on the right gives 

ZnFn+Zi2Fn= / 

(A4.3) 
Z21Y11 + ZasFn = 

Z21F12 -\- Z22Yi2 = ■* 

Transposing the matrices in these equations leads to other relations. Thus, 
from the first we obtain FuZn + Y12Z21 = I. These equations also yield 
expressions for the F's in terms of the Z's and vice versa. 

A somewhat similar treatment involving equations (2.1) and (2.3) leads 
to expressions for the Z's in terms qIA,B,C and D. The F's may be like- 
wise expressed. These relations are given in the following table. 

Fii = DK^^ Yn = Zn - ZuZI^Z^ 

F,, = C - DB'^A - -B'~' YT2 = Z21 - Zi^ZTiZn 



(A4.4) 



F21 = -B i'ai = '^n - yr.n-^21^22 

F22 = B^^A Y22 ^ ^22 — Z2iZn Z12 

Zn = AC-' ZTi - Yn - YUYT2Y21 

Z12 = AC-'D - B = C'-' ZTi = F21 - Y22YV2YU 

Z21 = C-' Z^l = Yn - FnFn'Fw 

Z22 = C-'D Z22 = F22 - FaiFTi'Fia 

A = Z11Z21 — — 1^1 F22 

B = ZhZ^iZqo — Z12 — — ITi 

C = Z^^ = Yu - YriYTlY22 

D — Z21 Z22 = — F11F21 

r^ BT' r ^' -^'1 ^^' -^C' = I AB' = BA' 

[c d] ~l-C A'] CD' = DC' DA'-CB' = I 
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The following equations in which X is an arbitrary scalar multiplier may 
be verified by equating coefficients of powers of X and using the relations 
just given. 

(Zai - XZu)(Fii + XF12) = (\Zi2 ~ Z22){\Y22 + Y21) 

{\^Zn - XZn - XZ22 + Zni)(Fn + XF12) 

= {\Zn - Zn){\n\2 + XFn + XF22 + F21) (A4.5) 



r-F2i irx^ -/ \B n r 

L Z21JL XC XZ)-/J L 



XF22 + F21 x/ 

X/ XZ22 — Z21 



Sometimes it is of interest to obtain the elements of F12 , say, when 
Zii , Z22 , Fn , F22 are known. Relations helpful in studying this problem 
are 

FiiZiiFi2 = Fi2Z2.F22, YuZnY^^- Yu= VuZnY^i 

YnY^tYn = r„ - ZTi Z12 = -ZiiFi2Fr2 

FaiFTi'Fi. = F22 - Z22 Z21 - -I'72'(FiiZn - /) 

When the typical section is symmetrical some simplification takes place 
and we have 



F,i = F.. 


Zu = Zo. 


A = D' 


AB = BA' 


F,2 = F21 


Zl2 = Z2I 


B = B' 


A'C = CA 






c = a 


A^~BC = I 



(A4.6) 

ZiiFu + ZuYn = I A'B'U ~ C = B'' 
ZnYu + Z12F,, = 

APPENDIX V 
Properties of the Matrix GX- + H\ + C 
The matrix 

/(X) = 6V + Z/X + G' (A5.1) 

which entered the discussion of the case of unsymmetrical sections has a 
number of interesting properties which are given below. G and H are 
square matrices with m rows each, and H is required to be symmetrical. 
As before, we shall denote by Xi , ■ ■ - Xm , X7\ ■ ■ ■ X^' the 2 m roots of the 
determinantal equation 

|/(X) 1 ^ 



.. .^ ,, -., . -^.«*„^^gt?«B;^— . . - '-^Lt~-^^'^aig.^^-^^^XJ. 
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and we shall suppose these roots to be distinct. Let the column kr and the 
row Ir be such that 

krlr = F{\r) (A5.2) 

where F{\) is the matrix adjoint to/(\), and let the square matrices K and 
L be defined by 

h 



K = [ki,k2, '--kj, L = 



U 



L.1^- 



CA5.3) 



Comparison of (A5.3) and (2.50) suggests that when G and H are expressed 
in terms of the F's we have the relations 

K = P, L=T' (A5.4) 

The method of choosing the column pr and the row pr shows that they are 
related by 

prVr = yrFCKr) 

instead of (A5.2) where 7r may turn out to be any non-zero constant, and 
consequently equations (A5.4) are not satisfied in genej^al. Nevertheless K 
and L may be regarded as particular choices for P and P'. In the same way 
K and L may be regarded as particular choices for Q and Q' when G and H 
are expressed in terms of the Z's. There is still some arbitrariness connected 
with K and L since the product kX is unchanged when the K is multiplied 
by some number and h is divided by the same number. 
The relations which correspond to (2 =^2) and (2.57) are 

GKA^+HKA+G'K = 

(AS. 5) 

GL'A"' + Z/i'A"' + G'L' = 

where A is the diagonal matrix whose e'ements are Xi , X2 , ■ • • Xm ■ These 
relations are consequences of the properties of kr and h. Two more rela- 
tions may be obtained by transposition. From the first of (A5.5) and the 
transposed of the second it follows that 

GKAKT' + ff + G'KK-'K-' =0 

(A5.6) 

L~'ALG + H + L~'A~'IG' = 

where it is assumed that the reciprocal matrices K~^ and L exist. Com- 
binations similar to KAir\ KA^'K'', etc. enter the expressions (2.59) 
for Zo and Yo. 
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By differentiating the equation 

/(X)F(X) = A(X)/, 
where A(X) is the determinant 

A(X) = |/(X)| = |GX'+ffX + G'|, 
it may be proved that 

GKAK"' + // + 7~'ALG = L~^EK~^ 
CKAT^FT' + /f + L~^K~^LG' - -L~'EK~^ 
in which E is the diagonal matrix 

r (1) 
A(Xi) ■ ■ ■ 
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E = 



(1) 
A(X2) 



(1) 
A(X™)J 



and 



w r 

A(X,) = [ 



(1) 



Since the roots X^ are assumed to be distinct, A{\r) ^ 0. 
We also have the equations 

KE~'L = L'ET^K' 
GKAE~'l - GL'A-'sr'K' - 7 



(A5.7) 



{A5.8) 



The first equation of (A.S.8) shows that KET^L is a symmetrical matrix. 
From this and the second equation it follows that 

GKAK~^ - GL'A~^L'~^ = L~'EK~^ 



(A5.9) 



From the first of equations (A5.7) and the second of {A5.6) 
GK\K~^ - L~'a~'LG' = L-'EK~' 



(A5.10) 

and the comparison witii {A5.9) shows that the matrix GL'A~ V~ is sym- 
metrical. The other matrices of this type are also symmetrical as may now 
be seen from equations (A5.6) and (A5.7). Results of this sort may be 
obtained from physical principles by noting that Zg and To must be sym- 
metrical matrices. 



,-4- -._^.,^.- 



^• ii il> i >i i*ti "iu. . \_. nu}tia»«fc;^- 
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As an example of the application of these formulas we assume that G and 
H are expressed in terms of the F's. Then we may take K and L' to be 
particular choices of P and P and equation (A5.9) becomes 

FisPAP"^' - Fi2pA-'p-' ^ (PT'Ep-\ 

From equations (2.59) and (2.62) 

To +Yo= Fi2^AP-' - F,2pA-'P-' 

and hence 

P'{Yo + Yo)P - E. 

For the more general choice of P and P allowed in §2.8 the diagonal matrix 
E is replaced by a general diagonal matrix. Similarly it follows that 

Q'(Z, + Z,)Q 

is a diagonal matrix. 



